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Abstract 

Classical and quantum mechanical results are reported for the single particle 
motion in a harmonic oscillator potential which is characterized by a quadrupole 
deformation and an additional octupole deformation. The chaotic character of the 
motion is strongly dependent on the quadrupole deformation in that for a prolate 
deformation virtually no chaos is discernible while for the oblate case the motion 
shows strong chaos when the octupole term is turned on. 
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The mean field approach is of central importance in any theoretical description of a 
many body system. In nuclear physics it is the basis for the shell model and its extensions 
such as collective states. Its success has been further demonstrated in the application to 
deformed nuclei and to metallic clusters where spherical symmetry is given up because of 
experimental evidence. Usually, quadrupole deformation is considered to be the major de- 
viation from spherical symmetry. However, more recently a possible octupole contribution 
has been taken into account for a number of reasons 0. 

Inclusion of an octupole term in addition to a quadrupole term renders the classical 
single particle motion nonintegrable. In fact, the system turns out to be chaotic. This has 
been discussed by a number of authors with various degrees of simplification P]. A nicely 
systematic approach is given in [^], where the motion in a quadrupole deformed cavity is 
analyzed and the terms that give rise to actual chaotic behaviour are clearly distinguished. 
In a recent investigation ||^ the study of classical motion in a cavity with oscillating walls 
of even and odd higher order multipoles has led to interesting conclusions about elastic 
versus dissipative behaviour of a noninteracting gas depending on the integrability or 
nonintegrability of the equations of motion. Surfaces of section are used to discern the 
onset and degree of chaotic motion. 

This paper is similar in spirit in that we investigate a simplified model where we leave 
out terms, which, albeit physically important, are prone to blur the analysis when the 
interest is focussed on the essentials that give rise to chaotic behaviour. Since our aim 
is directed not only to the classical but also to the corresponding quantum mechanical 

motion, we leave out the spin-orbit term and the / -term pr Nilsson model to render as 
closely as possible the analogy between the classical and quantum case. The importance 

of the I -term is well known in nuclear physics and for metallic clusters^. To determine 
its role in the context of chaotic motion the corresponding classical case ought to be 
studied. We defer treatment of this term. Despite the simplifications, such a model 
allows to understand the main features of shell structure effects, for example, in super- 
(hyper) deformed nuclei |0, §]. Recent experimental data of superdeformed /^-isomers in 
nuclei and electronic shell structure effects in metallic clusters clearly underline 
the importance of oblate deformation. Therefore we investigate in the present paper the 
effect of the octupole term for prolate and oblate deformation; the oblate/octupole case 
has not been dealt with in ref . . We analyse the case of zero temperature which is good 
for nuclei. For metallic clusters finite temperature should be considered |T^; however, our 
interest is focussed on shell structure whose character is unaffected by temperature except 
for the amplitude. Only axially symmetrical terms are taken into account which brings 
down to two the number of degrees of freedom of the motion. 
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The single particle motion is considered in the potential 

vi.s.A = —[i^^,^^^^) (1) 

where (? = ^ in cartesian coordinates x,y,z. We recognize an axially symmetric 
harmonic oscillator with frequencies Ux = ujy = bUz with an octupole term, written in 
the cylindrical coordinates {g,z,(f)). For the results presented the parameters are chosen 
such that the levels are 15MeV apart for 6 = 1 and A = which corresponds to a three 
dimensional isotropic harmonic oscillator. For A = and 6 > 1 (6 < 1) we have the mean 
field potential of a prolate (oblate) nucleus. Note that, choosing a suitable set of different 
parameters, we deal with a metallic cluster. The octupole term (A 7^ 0) is proportional 
to r'^Y^Q with 1^30 being the third order spherical harmonic and r'^ = g'^ + z"^. 

If |A| < Ac we are dealing with a proper bound state problem. Here Ac is defined to be 
the value for which the potential no longer binds, for |A| > Ac the potential tends to —00 
along one or two directions. The direction and the value of Ac depend on the quadrupole 
deformation b. For prolate nuclei (6 > 1), Ac = 1/(26^) and the potential opens its 
valley along the positive (negative) z-direction for negative (positive) A. For oblate nuclei 
(0.5 < b < 1) the other possible direction which is along the line g = sign(A)Q;2;, with 
a ~ 0.4, is of increasing importance. At the value b ~ 0.58 valleys along the two directions 
^ = and g = OAz open simultaneously for Ac ~ 1.5 while for a still smaller value of b, 
say for b = 0.5, the valley along the direction g = OAz opens for Ac ~ 1.64 while the one 
along the 2;-axis now opens for a larger value of A. (Analytic expressions for Ac and the 
direction a as functions of b exist but are of little interest.) Chaotic motion is expected 
to become the more pronounced the nearer the parameters for a bounded motion are 
to those for an unbounded motion. Thus increasing chaotic behaviour is expected when 
A approaches Ac. Furthermore, since also the geometry of the potential for A = Xc{b) 
depends on b, the chance of an escape for A > Ac also depends on b. This chance reflects 
upon the amount of chaos, i.e. the Lyapunov exponent, for the bounded motion prevailing 
for values A < Ac(&). In this way the amount of chaos is expected to be greatest for 
b ^ 0.58 where two valleys can serve as an escape route if A = Ac. In addition, at 6 = 2 
(superdeformed prolate nucleus) less chaos is expected than at 6 = 0.5 (superdeformed 
oblate nucleus), since the opening of the phase space happens through a small bottle neck 
in the former case thus considerably reducing the chance for escape, while the chance for 
escape at 6 = 0.5 is enhanced as the phase space opens more widely. 

The results of the numerical integration of the equations of motion confirm all the 
expectations. Axial symmetry yields the constant of motion = g'^(f) of the three dimen- 
sional motion. We present results only for = . A non-zero value of the z-component 
of the angular momentum does not produce new insights. Since the potential scales as 
V{'yr) = 7^V"(r) it suffices to investigate one energy only |]TT]. 



For the superdeformed prolate nucleus (6 = 2) there is hardly any chaotic behaviour 
discernible in the classical motion for all A < Ac. When looking at surfaces of section 
which we have taken at g = (recall = 0) in the (p^ — 2;)-plane there is some scattering 
of the dots when A is very close to Ac, but the Lyapunov exponent is virtually zero. There 
is, however, a proliferation of periodic orbits, an aspect important below when we discuss 
quantum mechanics. This is illustrated in Fig. (la) where surfaces of section are displayed 
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for 6 = 2 and A = 2/3 Ac. 

For decreasing values of b the motion becomes increasingly chaotic up to the maximally 
chaotic case at 6 ~ 0.58. Surfaces of section are illustrated in Fig. (lb) for A = l/3Ac. 
Note, that for b = 0.58, no structure would be discernible for A = 2/3Ac. We have 
compared results for A = pXc{b) with jS = 0.2, 0.4, 0.6 and 0.9. The trend is uniform 
in that the Lyapunov exponent shows the behaviour as indicated in Fig. (2). The figure 
presents values referring to chaotic orbits and avoids initial conditions within regions of 
stability. Such regions of stability gradually disappear when A approaches Ac if 6 < 1. 

We find the expected proliferation of periodic orbits. However, we refrain from dis- 
cussing them in great detail as the aspect of the possible retrieval fl^ of periodic orbits 
from the quantum mechanical spectrum is discussed in a forthcoming paper. 

The quantum mechanical treatment is straightforward in principle. In the spirit of 
previous work [1T3 we use for the full problem which is of the form Hq + XHi a repre- 



sentation where Hq is diagonal. The basis chosen is referred |0] to as the basis using the 
asymptotic quantum numbers n±, Uz and A where n± = n+ + n_. For a fixed value of A 
this leaves two quantum numbers (reflecting the two degrees of freedom) to enumerate the 
rows and columns of the matrix problem. For A = the diagonal entries of Hq are thus 
En^,n^ = huj{n_i_ + 1 + {uz + l/2)/b). The matrix elements of Hi are obtained from those 
of 2; ~ (aj -|- Qz) and ~ (A+a!,_ + + A^A^ + A\_AI.), where, in terms of the usual 

boson operators a^^y, A± = (a^ ^iay)/\/2. The effect of truncation was tested by looking 
at the variation of the lower end of the spectrum when the dimension of the matrices was 
increased. There is certainly a dependence on b and A. For A < 0.9Ac and 0.5 < b < 2 the 
variation was less than 1% for the first 300 levels obtained from 1600x1600 dimensional 
matrices. 

For demonstration we display parts of spectra in Figs. (3). It is obvious, and in fact 

quantitatively confirmed in Figs. (4) where the relevant statistical analyses are shown, that 

the quantum mechanical results are in line with the classical cases. The level repulsions 

in the superdeformed prolate case are very weak thus giving rise to a nearest neighbor 

distribution (NND) which appears nearer to an integrable case than to the typical Wigner 

distribution. Of physical interest are the pronounced new shell structures that emerge 

near to A = 0.5Ac and A = 0.65Ac in the superdeformed prolate case. This pattern is 

directly related to the periodic orbits indicated in Fig. (la) as will be discussed in detail 

^2 

in a forthcoming paper. Since we have left out terms like spin-orbit coupling and / we 
cannot claim that such structure where we find them. However, the essential point is the 
fact that such structures will always emerge. Since the situation is so close to integrability 
for all A < Ac, the level repulsions are always weak; as a consequence, shell structures 
are bound to emerge for some values of the parameters. We note, that the consequent 
periodicity in the spectrum will produce sharp lines in the Fourier transform of the level 
density fl^, a reflection of the many periodic orbits [1^ found in the corresponding 
classical case as was mentioned above. 

These findings nicely contrast with the oblate case where the spectrum and in partic- 
ular the NND have all signatures of chaos. For sufficiently large values of the octupole 
strength A all periodic structure is destroyed and there is no scope for new magic numbers. 
We mention that the NND best approaches the Wigner surmise for b ~ 0.58. 
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In summary, we have found for a problem typical for nuclear physics and for metallic 
clusters that chaotic behaviour should be expected in principle even for the single particle 
motion if deformations of higher order than quadrupole are taken into account. Moreover, 
for the prolate and in particular the superdeformed prolate case, there is a remarkable 
stability against chaos when octupole deformation is switched on. This result is in agree- 
ment with a prediction of an octupole instability of super- and hyperdeformed nuclei |TB 



based on realistic calculations. We conjecture that this pattern prevails also when axial 
symmetry is broken, i.e. when terms of the form Y^m, m ^ are taken into account. To 
what extent stability against chaos can be associated with stability in general terms, - it 
is a well known fact that there are more prolate than oblate nuclei -, is subject to further 
investigation. In the spirit of ref. we should expect prolate nuclei to behave rather 
elastically in contrast to oblate nuclei where more dissipative behaviour is anticipated. 
We may speculate that the absence of shell structure for the oblate/octupole case could 
prevent the existence of stable oblate/octupole deformed clusters. 
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Figure Captions 



Fig.l Surfaces of section at = for 6 = 2, A = 2/3Ac (left) and b = 0.58, A = l/3Ac 
(right). In the left part stability islands are clearly discernible for winding numbers 2:5, 
1:2 and 4:7. The right part is dominated by chaotic motion; some of the remaining islands 
are indicated by ellipses. 

Fig. 2 Lyapunov exponent as a function of A/Ac. The lower solid line is for b = 1, the 
upper solid line for 6 = 0.58 and the dashed line for b = 0.5. The values for 6 = 2 are too 
small to appear on the diagram. 

Fig. 3 A section of energy levels as a function of A/Ac for 6 = 2 (top) and b — 0.5 
(bottom) . 

Fig.4 Nearest neighbor distribution of the spectra shown in Fig. 2 for A/Ac ~ 0.5 as 
a function of the unfolded energy. For the top (6 = 2) 200 levels and for the bottom 
(6 = 0.5) 600 levels have been used. 



6 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/94060 1 2 v 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/94060 1 2 v 1 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/94060 1 2 v 1 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/94060 1 2 v 1 



